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Abstract 

In the present work, the bulk degradation of SiC in hadron (pion and pro- 
ton) fields, in the energy range between 100 MeV and 10 GeV, is characterised 
theoretically by means of the concentration of primary defects per unit flu- 
ence. The results are compared to the similar ones corresponding to diamond, 
silicon and GaAs. 
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1 Introduction. General properties 

Silicon carbide is a new interesting material for use in optoelectronic and elec- 
tronic devices, detectors, for the use in nuclear waste technology or in fusion re- 
actors. In all these applications, this semiconductor material could be exposed to 
extreme conditions of temperature, to chemical agents and to radiation fields. 

SiC is not a simple semiconductor, it is a whole class of semiconductors because 
of its polytypism [[I). Polytypism refers to one-dimensional polymorphism, i.e. the 
existence of different stackings of the basic structural elements along one direction. 
More than 200 polytypes have been reported 0. Only few of them have practical 
importance. These include the cubic form 3C(/3), and the AH and 6H hexagonal 
forms. For the cubic polytype, the symmetry group is Tj, while for the hexagonal 
ones this is C\ v . 

Some material characteristics of different polytypes of SiC, useful for its utilisa- 
tion in radiation fields, are presented in comparison with the corresponding ones of 
diamond, silicon and GaAs in Table 1. The compilation of data has been realised 
using Refs § | | | @, |, |, 0, 0, |§ . 

The listed material properties could explain why electronic devices made from 
SiC can operate at extremely high temperatures without suffering from intrinsic 
conduction effects due to its wide bandgap. Also, this characteristic allows to emit 
short wavelength light when fabricated into a light emitting diode structure. 
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The high breakdown electric field, approximately one order of magnitude higher 
than the silicon or GaAs ones, and the high thermal conductivity, enable the fabri- 
cation of very high-voltage, high-power devices that dissipate large amounts of the 
excess heat generated. Additionally, it allows the devices to be placed very close 
together, providing high device packing density for integrated circuits. SiC devices 
can also operate at high frequencies (microwave) because of the high saturation 
electron drift velocity in the material. 

Some SiC material characteristics are useful for detection applications: its ra- 
diation length lies between the corresponding values of diamond and silicon; its 
relatively wide band gap (2.2 3.3) eV and a relative dielectric constant of 7.04 al- 
lows the assembly of detectors with low noise. With low mobilities: 450 cm 2 V _1 s _1 
for electrons and 115 cm 2 V _1 s _1 for holes respectively, SiC does not permit a fast 
signal development and therefore is not able to operate in high fluxes of particles. 

The crystalline materials for semiconductor devices used in high fluences of par- 
ticles are strongly affected by the effects of radiation. Up to now, in spite of the ex- 
perimental and theoretical efforts, most of the problems related to the behaviour of 
semiconductor materials in radiation fields, the identification of the induced defects 
and their characterisation, as well as the explanation of the degradation mechanisms 
are still open problems. 

Diamond has the reputation of being a radiation hard material and is considered 
as a competitor to silicon. It is supposed that the concentration of primary defects 
induced in SiC in radiation fields is between the corresponding ones in diamond and 
silicon respectively. 



2 Mechanisms of degradation 

At the passage of the incident charged particle in the semiconductor material 
same of its energy is deposited into the target. The charged particles interact with 
both nuclei and electrons in a solid. The total rate of energy loss, could, in general, 
be divided artificially into two components, the nuclear and the electronic part. 

The energy lost due to interactions with the electrons of the target gives rise to 
material ionisation, while the energy lost in interactions with nuclei is the origin of 
defect creation. The primary radiation defects are vacancies and interstitials. 

There is not one physical quantity dedicated to the global characterisation of 
the effects of radiation in the semiconductor material. A possible one is the concen- 
tration of primary radiation induced defects per unit particle fluence (CPD) fl3|| . 
It is not proportional to the modifications of material parameters due to the irradi- 
ation because of the time evolution of vacancies and interstitials - see the discussion 
in reference ]l3. Due to the fact that it does not exist a theory describing these 
processes of interaction, phenomenological models are used. 

In this paper, the mechanism considered in the study of the interaction between 
the incoming particle and the solid, by which bulk defects are produced, is the 
following: the particle, heavier than the electron, with electrical charge or not, 
interacts with the electrons and with the nuclei of the crystalline lattice. The 
nuclear interaction produces bulk defects. As a result of the interaction, depending 
on the energy and on the nature of the incident particle, one or more light particles 
are produced, and usually one or more heavy recoil nuclei. These nuclei have charge 
and mass numbers lower or at least equal with those of the medium p5[ . After the 
interaction process, the recoil nucleus or nuclei, if they have sufficient energy, are 
displaced from the lattice positions into interstitials. Then, the primary knock- 
on nucleus, if its energy is large enough, can produce the displacement of a new 
nucleus, and the process could continue as a cascade, until the energy of the nucleus 
becomes lower than the threshold for atomic displacements. Because of the regular 
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nature of the crystalline lattice, the displacement energy is anisotropic. The primary 
interaction between the hadron and the nucleus of the lattice presents characteristics 
reflecting the peculiarities of the hadron, especially at relatively low energies. If the 
inelastic process is initiated by nucleons, the identity of the incoming projectile is 
lost, and the creation of secondary particles is associated with energy exchanges 
which are of the order of MeV or larger. For pion nucleus processes, the absorption, 
the process by which the pion disappears as a real particle, is also possible. 

The energy dependence of cross sections, for proton and pion interaction with 
the nucleus, present very different behaviours: the proton-nucleus cross sections 
fl6| , [It]] decrease with the increase of the projectile energy, have a minimum at rel- 
atively low energies, in the range of 100 -j- 200 MeV, followed by a smooth increase, 
while the pion nucleus cross sections present for all processes a large maximum, at 
about 160 MeV, reflecting the resonant structure of interaction (the A 33 resonance 
production), followed by other resonances, at higher energies, but with much less 
importance Jl8|]. Because of the multitude of inelastic channels, in the concrete eval- 
uation of defect production, the nuclear interactions has been modelled previously 
in some simplifying hypothesis, see refs. 13 lj, 19, po| . 



The process of partitioning the energy of the recoil nuclei (produced due the 
interaction of the incident particle with the nucleus, placed in its lattice site) by 
new interaction processes, between electrons (ionisation) and atomic motion (dis- 
placements) is considered in the frame of the Lindhard theory pl[ |. 

3 Results and discussions 

The concentration of the primary radiation induced defects per unit fluence in 
SiC has been calculated using the explicit formula: 



r Y,{%) L(E Rl)sicd n (i) 
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where E is the kinetic energy of the incident particle, Nqi^ is the atomic density of 
the C{Si) in SiC, Acfsi) is the atomic number of the C (Si), E^ciSi) - the threshold 
energy for displacements in the C (Si) sublattice of SiC, Em - the recoil energy of 
the residual nucleus produced in interaction i, L {Em) - the Lindhard factor that 
describes the partition of the recoil energy between ionisation and displacements 
and dai/dfl - the differential cross section of the interaction between the incident 
particle and the nucleus of the lattice for the process or mechanism i, responsible in 
defect production. The atomic density Nc(si) in SiC is a material parameter which 
depends to the polytypc through structural constants. The formula gives also the 
relation with the non ionising energy loss (NIEL). Na is the Avogadro number. It 
is important to observe that there exists a proportionality between the CPD and 
NIEL only for monoelement materials, and this proportionality can be extended 
only approximately for compounds with remote elements in the periodic table, e.g. 
for GaAs. 

The CPD allows the comparison of the effects produced by the same particle in 
different materials, while NIEL is especially used for the comparison of the effects 
produced in the same material by different particles. 
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The partition of the energy of the incoming particle in the semiconductor mate- 
rial between ionisation and displacements has been calculated analytically, solving 
the general equations of the Lindhard theory in some physical approximations. De- 
tails about the hypothesis used could be found in reference |^2|. In the particular 
case of the SiC, a binary compound, the Lindhard curves have been calculated sep- 
arately for each component, and the average weight Bragg additivity has been used. 
The energy spent into displacements is represented versus recoil energy for different 
recoils, characterised by their charge and mass numbers, in SiC, in Figure 1. 




Figure 1: The displacement energy versus the energy of recoils (characterised by their 
atomic and mass number) in SiC. 

All curves start, at low energies, from the same curve; they have at low energies 
identical values of the energy spent into displacements, independent on the charge 
and mass number of the recoil. At higher energies, the curves start to detach from 
this main branch. This happens at lower energies if their charge and mass numbers 
are smaller. The maximum energy transferred into displacements corresponds to re- 
coils of maximum possible charge and mass numbers (corresponding to the heaviest 
element, Si in this case). The curves present then a smooth increase with the en- 
ergy. For the energy range considered here, the asymptotic limit of the displacement 
energy is not reached. 

For the cross sections for proton and pion interactions with carbon and silicon, 
experimental data and theoretical extrapolations have been used - see references 
JTg| , |23| , [24| and references cited therein. The cross sections for pion Ga and As 
nuclei are not available from experiment. Due to this fact, parametrisations of 
the energy dependence of the cross sections have been tried. Details about the 
procedure are given in reference |25| . 

As stated in Section 2, the minimum energy required to form a stable vacancy - 
interstitial pair is usually called the threshold energy for displacements (Ed). The 
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regular nature of the crystalline lattice conduces to its anisotropy, atoms being more 
readily displaced along some directions than in others. This anisotropy is also the 
origin of differences in the rate production of defects as a function of orientation. In 
the literature, both calculated and experimental values for are reported for different 
crystals. 

In principle, for SiC, the situation is even more complicated and a threshold for 
atomic displacements could be defined for each polytype, for each atomic species, 
and for each direction in the crystal. 

Some calculations have been reported in the literature for Ed', in SiC, along some 
directions of interest, by classical molecular-dynamics simulations using empirical 
potentials |2(| , and also using modern first-principle methods |27]] . 

Experimental values are obtained from high voltage electron microscopy |2q| , 
energetic ion beam spectroscop y |pS |, and Rutherford backscattering analysis along 
different channelling directions |lC| , 29[ [j0[ . 

In the opinion of some authors J30J , all these methods give consistent results, the 
recommended mean threshold energy for atomic displacements being of the order 
of 35 eV on the Si sublattice, and of the order of 20 eV on the C sublattice. These 
are the values we used in the calculations. 

The energy dependence of the CPD calculated (as stated above) for pions in 
SiC, together with the estimated values for proton irradiation is presented in Figure 
2. In the same figure, for comparison, the CPD produced by pions and protons in 
diamond |lj| |24j], silicon and GaAs are also represented. The values corresponding 
to protons in silicon are averages from the calculations reported in references ]l2| 
and pl| |. The curve for protons in GaAs is obtained from the results reported in 
reference and is estimated by extrapolation in the energy range of interest for 
the present study, while the curve for pions in GaAs is from reference [^3| . It could 
be observed that the degradations induced by pions and protons are very different 
and reflect the peculiarities of the interactions of the two particles with the nuclei 
of the semiconductor. For pions, the maximum of degradation is produced in the 
region of the A33 resonance, corresponding to about 140 160 MeV kinetic energy. In 
comparison with this behaviour, the maximum of the degradation is produced, for 
protons, at low energies, and a smooth and slow decrease of the CPD is calculated 
with the increase of the energy. For the energy range considered in this paper, 
with the exception of the resonance region, the degradation induced by protons is 
higher than that for pions in SiC, but for both incident particles, the concentration 
of primary defects produced per unit fluence are closed between 30 and 40 cm . 
The degradation induced by hadrons in SiC is between the diamond and silicon 
degradations. 

Also, it could be observed that the pion degradation becomes more important 
with the increase of the mass number of the material, in the whole range of energy 
investigated. 

The CPD for 3C, 4H and 6H polytypes of SiC has been evaluated separately, 
considering the structural characteristics of each polytype, and the same displace- 
ment energies for all of them; no significant differences between their behaviour in 
radiation fields has been found. These theoretical results are in agreement with the 
experimental data reported in reference |32]. 



4 Summary 

In the present paper the concentration of primary induced defects per unit par- 
ticle fluence has been calculated for the case of pion and proton irradiation, and 
the effects of these particles have been found intermediate in SiC between the cor- 
responding ones in Si and diamond. 
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Figure 2: Concentration of primary induced defects on particle unit fluence for pions 
(continuous) and protons dashed) in (up to down): gallium arsenide, silicon, silicon carbide 
and diamond. 



Negligible differences have been found between the 3C, 4H and 6H polytypcs 
in what regards the effects of pion and proton irradiation. 

The Lindhard curves characterising the partition of the recoils energy between 
ionisation and displacements have been calculated using analytical approximations. 

The main properties of SiC useful for its utilisation in radiation fields have been 
reviewed. 
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